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Superposition of gravitational waves generated by astrophysical sources is expected to give rise 
to the stochastic gravitational-wave background. We focus on the background generated by the 
ring-down of black holes produced in the stellar core collapse events across the universe. We sys¬ 
tematically study the parameter space in this model, including the most recent information about 
the star formation rate and about the population of black holes as a function of redshift and of 
metallicity. We investigate the accessibility of this gravitational wave background to the upcoming 
gravitational-wave detectors, such as Advanced LIGO and Einstein Telescope. 
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INTRODUCTION 


in Europe [61 1. 


The stochastic gravitational-wave background 
(SGWB) is expected to arise from the superposition 
of gravitational waves (GWs) from many uncorrelated 
and unresolved sources. Numerous cosmological and 
astrophysical models have been proposed. Cosmological 
models include inflationary models EHl , models based 
on cosmic (superWrings j5l4lCll|. and models of alterna¬ 
tive cosmologies [ll|, Astrophysical models (see 131 
for a review) integrate contributions from astrophysical 
objects across the universe including compact binary 
coalescences (CBC) of binary neutron stars (BNS) or 
binary black holes (BBH) Il4l424l. rotating neutron stars 
(NSs) magnetar s l34l44cl| . the first stars [dlj, 

and white dwarf binaries [42l| . 

Several searches for the isotropic ^ | and 
anisotropic SGWB [il, have been conducted us¬ 
ing data acquired by the first generation interferomet¬ 
ric gravitational-wave detectors LIGO 5^, 5l| and Virgo 
(5^ . These searches have established upper limits on the 
energy density in the SGWB, and have started to con¬ 
strain some of the proposed models [s, [2S[5i,[53. The 
second generation of gravitational-wave detectors is cur¬ 
rently beii^ commissioned, including Advanced LIGO 
(aLIGO) [5I, Advanced Virgo [H^, GEO-HF [H^, 
and KAGRA [H, |6^ . These detectors are expected to 
produce first data in 2015, and their strain sensitivity is 
expected to be 10 x better relative to the first genera¬ 
tion detectors. The third generation gravitational-wave 
detectors are also being conceptualized, such as the Ein¬ 
stein Telescope for which the design study was completed 
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One of the promising sources of gravitational waves, 
potentially detectable by the second and third genera¬ 
tion detectors, is the stellar core collapse process. The 
physics of core collapse is complex and it is expected to 

g iduce gravitational waves via several mechanisms 62- 
. While the full three-dimensional simulations that 
include all relevant processes are yet to be made, predic¬ 
tions have already been made about the gravitational- 
wave signals emitted during core collapse 62l-l64l| . These 
predictions have then been used to make estimates of the 
corresponding SGWB due to both standard and early 
(population III) stars [HI Ibbl - fT^ . These estimates are 
necessarily only approximate, both because it is currently 
not well understood how the gravitational-wave signal 
depends on the progenitor stellar parameters such as the 
mass or spin, and because the rate of core collapse events 
is uncertain. 

In this paper, we focus on the ringdown of the black 
hole (BH) following the core collapse. The GW spec¬ 
trum emitted by this process is relatively well under¬ 
stood, since most of the energy is dissipated via the ring- 
down of the I = 2 dominant quasi-normal mode l76l ) , 

as is also confirmed in simulations [silTl. The SGWB 
due to this source was first examined by Ferrari et al [65| . 
We revisit this model of SGWB, taking into account the 
most recent information about the star formation rate, 
about the population of black holes, and about the ef¬ 
fects due to varying metallicity. We also investigate the 
accessibility of the model to the upcoming second and 
third generation GW detectors. In Section 2 we present 
the general aspects of the calculation of the astrophysical 
SGWB, and we discuss the star formation rate used in 
this study. In Section 3 we discuss the SGWB due to 
the core collapse to black holes (Model 1)—this has been 
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discussed in the literature but we revisit it here using 
the latest star formation rate. In Section 4, we discuss 
Model 2, which builds on the Model 1 by including the 
effects of metallicity. In Section 5, we discuss Model 3 
which is based on a Monte Carlo simulation of the SGWB 
due to core collapse to black holes, taking into account 
the Star Track m numerical simulation of the black hole 
population. We summarize our results in Section 6. 


ASTROPHYSICAL GRAVITATIONAL-WAVE 
BACKGROUND 


To compute the SGWB due to astrophysical GW 
sources, we follow the formalism used in [isj, l20|, l37|. In 
particular, we define the normalized energy density in 
gravitational waves: 


^^Gw(/) 


I dpGW 

Pc d In / ’ 


( 1 ) 


where PGw(/) is the energy density in gravitational 
waves at the observed frequency /, and pc = g^°Q is 
the critical energy density needed to close the universe. 
Here, Hq is the present value of the Hubble parameter, 
taken to be 68 km/s/Mpc, G is the Newton’s constant, 
and c is the speed of light. The energy density can then 
be rewritten in terms of the integrated flux F(f) over 
redshift z: 


^Gw{f) = —F{f) 

PcC 

_ J_ f , Rzjz) dEgw 
PccJ ^47rr2(2;) dfe 


( 2 ) 


Here, r{z) is the proper distance, and dEcpj^i/df is the 
energy spectrum emitted by a single astrophysical source 
as a function of the emitted frequency /e in the source’s 
frame, /e = /(I + z). We will consider three different 
models for dEcy^/dfe below. The rate of astrophysical 
sources as a function of redshift, Rz{z), can be written in 
terms of the rate of sources per comoving volume Rv{z): 


Rz{z) 


Rv{z)^ 

dz 


Rv{z 


47rc 




F[q m, z) , 


(3) 


where E{Qm,^A, z) captures the dependence of the co¬ 
moving volume on redshift: 

E{VtMz) = \/Hm(I + z)^ Ha, (4) 


with Qm = 0.3 and ft a = 0.7 corresponding to the energy 
density in matter and dark energy respectively. Finally, 
the rate of astrophysical sources can be related to the 
star formation rate (SFR) 

Rv{z) = Abh 1 ^ ! (5) 

1 + z 


where Xbh captures the mass fraction of matter that 
ends up in progenitors of black holes. Bringing the above 
together, we get 

o (f\ _ SttG/Abh R*{z) dEcw 

J ‘^^^i + z)E{nM,nA,z) du ■ 

( 6 ) 


The star formation rate and its dependence on redshift 
have been studied extensively in the literature, and mul¬ 
tiple functional forms have been proposed 791484 1. His¬ 
torically, most of the star formation rate estimates have 
been based on luminosity measurements. We use the 
SFR derived in Behroozi et al. and Oesch et al. 
Hill. While luminosity measurements at redshifts up 
to ^ 2 are relatively well understood, measurements at 
high redshifts (up to z ^ 11) are subject to uncertainties 
due to the extinction by dust and due to the fact that 
early star formation takes place in faint galaxies that may 
be missed in magnitude-limited surveys. An alternative 
to luminosity measurements at high redshift is to use the 
rate of Gamma Ray Bursts (GRBs) 8^9^, which typi¬ 
cally results in a slower fall-off of the star formation rate 
at high redshifts. We use the GRB rate of [93 based on 
the normalization described in 


91|,|92|. 


As argued in [93| , the choice of the star formation rate 
has direct implications for the chemical and reionization 
history of the universe. The analysis presented in 
considered both the GRB and the luminosity data and 
arrived at two models of star formation rate that are con¬ 
sistent with the available metallicity data as well as with 
the reionization redshift and the optical depth measure¬ 
ment by WMAP 9^. Both models use the Springel & 


Hernquist functional form 


A* ( 2 ) = zz 


pe 


q{z-Zm) 


p 


-q + qePG-^ 


(7) 


with the following parameters [93j: 

• GRB-based model 93“ 9^: v = 

0.146 M0/yr/Mpc^ z™ = 1.72, p = 2.80, 
and q = 2.46. 

sH - 871: cont aining 

described hy u = 
0.178 MQ/yr/Mpc^, Zm = 2.00, p = 2.37, and 
q = 1.80, and population HI stars described by 
V = 0.00218 MQ/yr/Mpc^, Zm = 11.87, p = 13.81, 
and q = 13.36. 


• Luminosity-based model 
the normal mode stars 


Furthermore, both models use the same Salpeter initial 
mass function (IMF) (pirn) = where TV is a 

normalization constant defined such that 


(j){m)mdm = 1. 


( 8 ) 


Note that the normalization constant depends very 
weakly on the upper limit of this integral, so we set 
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FIG. 1: Comparison of the GRB-based (solid) and luminosity- 
based (dashed) models of star formation rate. 


TO 2 = oo. For the normal-mode stars we choose the 
lower limit of the integral to be mi = O.IMq. Since pop¬ 
ulation III stars are expected to be heavier, we choose 
mi = 36M0 for this population [o^. Figure [T] compares 
these two models of star formation rate as a function of 
redshift. We will compare our results for these two mod¬ 
els of star formation rate. Note that these two models 
predict substantially less star formation at high redshifts 
as compared with some of the models used in the past to 
estimate the stochastic gravitational-wave backgrounds 
(e.g. dH). 


MODEL 1 


We first consider the model of SGWB due to core col¬ 
lapse to black holes that was first proposed by Ferrari et 
al [g^. Numerical simulations 6^, 771 have shown that 
most of the energy is dissipated via the ringdown of the 
I = 2 dominant quasi-normal mode, whose frequency is 
given by 13, lii: 


i/*(m, a, a) 


A{a) 


am 

^(1-0.63(1-a)°-3) 


(9) 

( 10 ) 


where the mass of the black hole M is assumed to be a 
fraction a of the mass of the progenitor m (M = am), 
and a is the dimensionless spin factor ranging from 0 for 
a Schwarzschild black hole to 1 in the extreme Kerr limit. 
The energy spectrum of the single source is therefore 


= eamc^(5(/e - i/*(m,a, a)), (11) 

dfe 


where e is the efficiency of GW production, acting as a 
scaling parameter of the GW energy spectrum. Hence, 

^Gwif) = f dz dm (12) 

*5.ng J Jm„i„ 

R*{z) (j){m) m i5(/(l + z) — v^{m, a, a)) 

(1 -b z)E[^m, z) 


We can use the Dirac delta function to evaluate the in¬ 
tegral over redshift, resulting in; 


^Gw{f) 


SnCea 



Rt{z') (j){m) m 

“'™(l + 20i5(flM,f^A,20’ 

(13) 


where 


z = 


i, a, a) 


- 1 . 


(14) 


Note that since z' > 0 necessarily, the integral over m 
includes only the values of m for which i^*(m, a, a) > f. 

The range of integration in Eq. [T3] is defined by some 
minimum and maximum stellar masses that are expected 
to act as black hole progenitors. The lower end of this 
range, mmin affects only the high-frequency end of the 
GW spectrum, above ^ 200 Hz. Since most of the sensi¬ 
tivity of GW detectors to this SGWB comes from the low- 
frequency end of the spectrum, mmin does not strongly 
affect the accessibility of the model to the GW detectors, 
and we will fix it to 40Mq. Also note that am must 
be larger than ~ 4Mq, otherwise the black hole will not 
be produced. We verified that the GW spectrum is not 
very sensitive to the cutoff on am, for example setting 
the cutoff at 3Mq has negligible impact on Hqw- 

The high end of the integration range mmax affects 
the low-frequency end of the spectrum and could there¬ 
fore have a significant impact on the detectability of 
the model by GW detectors. Furthermore, there is 
currently much uncertainty in the largest mass of the 
black hole progenitor stars—for example, stars as mas¬ 
sive as ~ 200 — 3 OOM 0 stars we recently reported in the 
R136 star cluster in the LMC [9^. To capture this un¬ 
certainty, we will repeat our analysis for two values of 
Wmax = IOOM 0 and 5 OOM 0 to illustrate the importance 
of this high-mass cutoff. 

The free parameters of the model are therefore e, a, and 
a. Figure [ 2 ] shows example spectra for several choices of 
parameter values, in comparison with the expected sen¬ 
sitivity of Advanced LIGO and Einstein Telescope detec¬ 
tors. The e parameter is simply a scaling factor, and a 
effectively shifts the GW spectrum in frequency: low val¬ 
ues of a shift the spectrum to lower frequencies, at which 
the GW detectors are more sensitive. The parameter a 
has a more complex impact on the GW spectrum: low 
values of a shift the spectrum to higher frequencies and 
reduce its amplitude, both of which reduce the accessibil¬ 
ity of the model to GW detectors. Finally, Figure [5] shows 
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FIG. 2: The expected sensitivities of Advanced LIGO [9^ 
and Einstein Telescope [gJI detectors (assuming 1 year expo¬ 
sure) are shown in comparison with several examples of GW 
spectra obtained using Model 1 with e = 10~® and with the 
GRB-based model of star formation rate. While we assume 
co-located Advanced LIGO detectors in this study, their sen¬ 
sitivity to SGWB will be very similar to the sensitivity of the 
complete second-generation GW detector network, including 
Advanced LIGO, Advanced Virgo, GEO-HF, and KAGRA. 


the impact of the upper cutoff on the black hole progen¬ 
itor mass, rumax^ increasing this cutoff from IOOMq to 
5 OOM 0 extends the GW spectrum to ~ 5x lower fre¬ 
quencies, as expected based on Equations [TUI and [TTl 
We scan this parameter space, restricting 0 < a < 1 
and 0.1 < a < 0.5 (motivated by simulations [o^, as dis¬ 
cussed further below). For each point in this parameter 
space, we compute the spectrum flcwif) and compare 
it to the sensitivities of the Advanced LIGO and ET de¬ 
tectors, computing the likelihood function: 

L oc n (15) 


mation rates are shown. The curves in the e — a plane 
are decreasing with a, which is a consequence of the fact 
that increasing a pushes the GW spectrum to lower fre¬ 
quencies (c.f. Figure [2]), hence making the model more 
detectable by the GW detectors. Similarly, increasing the 
value of a pushes the GW spectrum to higher frequencies 
(c.f. Figure making the model less accessible to GW 
detectors and causing the increasing trend (with a) of the 
curves in the e — a plane. The upper cutoff on the black 
hole progenitors mass also has a significant impact on the 
accessibility of these models. As seen by comparing the 
two columns of Figure [3l changing rumax from IOOMq to 
500A/q lowers the sensitivity curves in the e parameter 
by 2-3 orders of magnitude. We verified that increasing 
nT-max to IOOOA/q further lowers the sensitivity curves by 
another factor of 2. 

The two star formation rates yield nearly identical pre¬ 
dictions, which is the consequence of the fact that the 
dominant contribution to LIgw comes from redshifts re¬ 
gion of 1-2, in which the two models of star formation 
rate agree well. It is also evident that the Einstein Tele¬ 
scope will provide a substantially better probe of this 
model than the second generation detectors. 

Finally, we note that the expected value of e in this 
mechanism for GW production is uncertain. While some 
of the past literature considers e as high as 0.01 [^ . 
previous simulations by Stark and Piran [771 gave an 
upper limit of e ^ 7 x 10“^ for an axisymmetric col¬ 
lapse. Accounting for more realistic scenarios, in par¬ 
ticular the pressure reduction that triggers the collapse, 
leads to e ^ 10“^ —10“® [o^. These values are clearly out 
of reach of the second-generation detectors as shown in 
Figure [21 but may be within reach of Einstein Telescope, 
especially in the case if massive black hole progenitors 
exist (corresponding to the case = 5 OOM 0 ). 


MODEL 2 


where the index i runs over frequency bins, Yi is the ex¬ 
pected measurement of the GW energy density in the 
bin i, cTi is the corresponding measurement error, and 
f^GW,i(ej O', ct) is the value of the energy density in the bin 
i for the given free parameters e, a, and a. For project¬ 
ing the future experimental sensitivities we set Y) = 0. In 
order to determine the accessibility of this 3-dimensional 
parameter space to future detectors, we marginalize (in¬ 
tegrate) the likelihood function over one of the parame¬ 
ters, and then compute the 95% confidence level contours 
in the plane of the remaining two parameters. 

Figure [3] shows the contours computed in this way for 
the € — a and e — a planes (top and bottom rows re¬ 
spectively) and for two choices of Wmax = IOOM 0 and 
5 OOM 0 (left and right columns respectively). Gurves for 
both Advanced LIGO and Einstein Telescope, and for 
both the GRB-based and the luminosity-based star for- 


We now modify the model discussed above to take into 
account the effects of metallicity Z. Metallicity has been 
shown to impact the maximum masses of the black hole 
progenitors as well as the fraction of t he p rogenitor mass 
that remains in the black hole ajZ) [97|. In particular, 
simulation results described in [97j | provide an explicit 
formula for the maximum black hole mass as a function 
of metallicity (Eq. 11 of 97|). Furthermore, this study 
shows the dependence of the remnant mass Mbh on the 
progenitor mass m for several values of metallicity - we 
model this dependence as linear, Mbh = a{Z)m, ex¬ 
tracting the value of a for several metallicity values, and 
then interpolating as necessary. For the sake of simplic¬ 
ity, we ignore the dependence of a on the progenitor mass 
m. The resulting curve for the average a{Z) is shown in 
Figure 01 Finally, we note that the minimum mass of 
black hole progenitors also depends on the metallicity, 
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m =100M 



m = 500 M 



m =100M 


m = 500 M 




FIG. 3: Top left: 95% confidence sensitivity contours of the Advanced LIGO and Einstein Telescope and for two models of star 
formation rate are shown in the e — a plane for Model 1, after marginalizing over the parameter a, and assuming mmax = IOOMq . 
Top right: Same as top left but for mmax = 500Mq. Bottom left: 95% confidence sensitivity contours of the Advanced LIGO 
and Einstein Telescope and for two models of star formation rate are shown in the e — a plane for Model 1, after marginalizing 
over the parameter a, and assuming mmax = IOOMq. Bottom right: Same as bottom left but for mmax = hOOM©. Note the 
significant change of vertical scale from models with mmax = IOOMq to models with mmax = 5OOM0. 


as we will see in Model 3—for Model 2, however, we 
will assume that mmin{Z) = 20 Mq, independent of the 
metallicity. Also, as in Model I, we require am > AMq 
in order for the black hole to be produced. 


fore rewrite the equation for 


^Gw{f) 


SnGfe 

~3Hf 


j dZip{Z-, z)a{Z) 


dz 


_ R*{z) _ 

(1 + z)E{^Im, z) 


(16) 


Metallicity is a function of redshift Z{z), and can be 
estimated using the chemical evolution code for the cho¬ 
sen star formation rate 93|. Figure [5] shows the mean 
metallicity evolution for the two star formation rates con¬ 
sidered in this paper. For any given redshift, we assume 
that the metallicity values range within a factor of 3 from 
the mean metallicity, with a uniform (flat) distribution 
denoted by '0(2'; z), and we average over it. We can there- 


^77Tmax(^) 


f TTlrnin (•^) 


dm 4>(m) m (5(/(l -I- z) 


l'^,{m, a, a)). 


We first solve the integral over the progenitor mass using 
the Dirac delta function: 




mo 


S(m — mo)a{Z)mm 
A(a) 

a(q) 

f{l + z)a{Zy 


^17) 

(18) 
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FIG. 4: Dependence of a (the fraction of the progenitor mass 
that remains in the black hole) on metallicity, extracted from 
the simulation results of [^ . 


Then, 


^Gwif) = 


SirCfe 


dz 


R*{z) 


dZ'ip{Z; z)a{Z)NmQ 


(1 + z)E(Q,m, z) 

0.65 


A(c 


0(f - /min)0(/max " /) (20) 

87rG/°-35e r R,{z) 

3il3A(a)0.35 J ^ (1 + z)i-e5E{nM, 

J dZi;{Z; z)a^-^HZ)NQ{f - /nnn)0(/max - /) 


If the condition in Eq. [TO] is not satisfied, the m-integral 
is simply zero. Note that N is the normalization of the 
initial mass function. 

Finally, if we ignore the metallicity dependence, the 
equation simplifies to: 



FIG. 5: Evolution of metallicity is shown for the two star 
formation rate models considered in this paper. Note that 

^Sun = 0.02. 


Clearly, the following condition must be satisfied for a 
non-zero signal: < mo < m^ax{Z). This trans¬ 

lates into 


_^(o)_ r __ 

TOmin(-^)(l + z)a{Z) m^^^{Z){l + z)a{Z) 

/max(^,a, z) > / > f^in{Z,a,z) (19) 


where in the last line we have defined /max and /min- 


^Gw{f) 


8TTGf -^^ea^-^^N f R^z) 

3 iJ3A(a)0-35 y (1 + zy-^^E{nM, ftA, z) 

0 (/ - fmin{a, z))0(/max(a, z) - /), (21) 


which is identical to Equation [TO] that was obtained by 
evaluating the z-integral instead of the m-integral. 


We note that the chemical evolution model that was 
used to determine the star formation rate 931 uses a value 


of a which varies between 0.08 a nd 0. 30 depending on the 
progenitor mass and metallicity |l00l |. However, we have 
confirmed that using the larger value of a = 0.5 in the 
chemical evolution model does not have a significant im¬ 
pact on the best ht for the star formation rate. In partic¬ 
ular, increasing a reduces the amount of metals released 
in the core collapse events by 10-20%, but the effect is 
well within the uncertainty on metallicity measurements. 


Hence, Model 2 is characterized by 2 parameters, e and 
a, which have similar effects on the GW energy spectrum 
as in Model 1. The third free parameter in Model 1, a, 
is now modelled as a function of metallicity and redshift 
based on the population synthesis simulations. Figure |6| 
shows a Model 2 spectrum in comparison with Model I 
spectra obtained for the same values of e and a and for 
the largest and smallest values of a. Model 2 spectrum 
can be seen as an effective average of the Model 1 spectra. 


Figure |7| shows the 95% confidence contours for the 
Model 2 in the e — a plane for the Advanced LIGO and 
Einstein Telescope detectors. Similarly to Model I, the 
two different star formation models yield nearly identical 
results, and the expected FT sensitivity is about 100 x 
better than that for Advanced LIGO. The contours are 
also about 2x lower in e than for Model I, which is a con¬ 
sequence of the fact that a{Z) is distributed differently 
in Models I and 2. 
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FIG. 6: Comparison of Models 1, 2, and 3, for e = 10 ® and 
a = 0.2. 



FIG. 7: 95% confidence sensitivity contours of the Advanced 
LIGO and Einstein Telescope and for two models of star for¬ 
mation rate are shown in the e — a plane for the Model 2. 


MODEL 3 


We now treat the dependence of a, rrimin, and rrimax 
on metallici ty more carefully, by using the results of 
StarTrack 101 . Il02l| . a sophisticated population syn¬ 
thesis code able to generate realistic populations of sin¬ 
gle and binary compact objects (neutron stars and black 
holes). The code is based on reyised formulas from |l03l | ; 
updated with new wind mass loss prescriptions, cali¬ 
brated tidal interactions, physical estimation of donor’s 
binding energy (A) and conyection driyen, neutrino en¬ 
hanced supernov a eng ines. A full description of these up¬ 
dates is given in |104l | . The two most recent updates take 
into account measurement s of initial parameter distribu¬ 
tions for massive O stars 105] as well as a correction of 
a technical bug that has limited the formation of BH-BH 


FIG. 8: Final mass of a remnant as a function of its pro¬ 
genitor’s mass and metallicity for single stellar evolution as 
predicted by our population synthesis model. Note that the 
maximum BH mass increases rapidly with decreasing metal¬ 
licity. Maximum BH mass for solar metallicity is only 15Mq, 
while for 1% solar metallicity it is SIMq. Our model also 
reproduces the observed mass gap between NSs and BHs. In 
our predictions there is no compact objects in mass range 
1.8 - 5.5Mq (Zq), 1.9 - 5.6Mo (0.1 Zq), and 1.8 - 6.7Mq 
(0.01 Zq). 


binaries for high metallicity (e.g., Z = 0.02). 

We evolve single stars until the formation of a compact 
object. Simulations are done for a dense grid of stellar 
metallicities: Z = 0.0001-0.03 with step of AZ = 0.0001 
(solar composition is Z = Zq rs 0.02). Two major fac¬ 
tors shape the initial (Zero Age Main Sequence) - final 
(compact object) mass relation; wind mass loss and core 
collapse/supernova compact object fo rmat ion. For wind 
mass loss we use O/B type winds from [l06] and for other 
evolutionary stages (e.g., LBV winds) formulae as cali- 
brat ed in [97|. We adopt the set of models presented by 


107] with the rapid core collapse/supernova mechanism. 


The explosion occurs within the first 0.1 — 0.2s driven 
by a convection and neutrino enhanced engine. This en¬ 
gine reproduce s 11081 the mass gap observed in Galactic 


X-ray binaries [1091 IllO] . The typical initial-final mass 
relations are presented in Figure [H] 

Pair-instability supernovae (PISNe) may disrupt the 
most massive stars without leaving behind a compact ob¬ 
ject. The original simulations indicate that stars in the 
initial mass ra nge M y,a,ms = 150 — 300Mq may be sub¬ 
ject to PISNe [111] (here, the subscript ”zams” stands 
for the ”zero-age mass sequence”). These initial esti¬ 
mates were translated into the final CO core mass (as 
this can be used for various metallicities and wind mass 
loss rates). The range for PISNe to disrupt stars with¬ 
out black hole form ation is found approximately in range 
Mco = 60— 130Mq [ 112 ] . For our assumed extent of IMF 































(up to 150 Mq) and the range of employed metallicities 
we are not in the regime of pair-instability supernovae 
(PISNe). However, recent discovery of stars as massive 
as ~ 200 — 3 OOM 0 in the R136 star cluster in LMC 


may possibly indicate that IMF extends beyond PISNe 
regime, implying that stars more massive th an 3 00Mf7^ 
may exist and form massive BHs (> IOOMq) [ll3l| . 

The evolutionary code StarTrack provides a(Z, m) for 
a grid of metallicity in the interval 0.0001 — 0.03, and pro¬ 
genitor mass between 7 — 150 Mq. It also gives mmin(.^) 
and minax(-^)) the minimal and the maximal masses to 
form a black hole for a given metallicity. Using these 
data, we can simulate a population of massive stars that 
undergo core-collapse to black holes. Our Monte-Carlo 
procedure is described below. 

We fix the spin parameter a and the efhciency e, and 
we proceeded as follows for TVmc = 10 ® sources : 


• We draw the redshift z from a probability distribu¬ 
tion constructed by normalizing the rate Rz{z) in 
the interval 0 — 20. 


Pz{z) 


Rzjz) 

N 


( 22 ) 


where N = Rz(z)dz and calculate the luminos¬ 
ity distance d^iz) = (1 + z)r(z). 

• We calculate the average metallicity Z at redshift 
2 and draw the metallicity from a uniform distri¬ 
bution in the interval \Z jZ — 3Z]. 


• By interpolating in the StarTrack data, we calcu¬ 
late mjain(Z) and mmax(^), as well as the corre¬ 
sponding mass fraction Xbh{Z) (c.f. Eq. [S|). For 
all metallicities, the maximal mass is assumed to be 
150 Mq while the minimal mass increases from 20 
Mq for Z = 0.0001 to 36 Mq for Z = 0.03. We also 
draw the mass of the progenitor from a distribution 
constructed from the initial mass function: 


Pm{m,Z) 




r^max(^) 
J 7Timin(^) 


4>(m)dm 


(23) 


We used two different models: the Saltpeter IMF 
introduced in the previous sections and a three 
component broken power-law IMF with slope of 
— 1.3 for initial mass Mzams = 0.08 — 0.5Mq, —2.2 
for Mzams = 0.5 - IMq, and -2.7 for M^ams = 

1 - 150M0 [TT^ . 

• We deduce a{Z,m) and then the mass of the BH 
Z) = a{Z,m)m, by interpolating in the 
StarTrack data. 


• We calculate the emission frequency 
iz*(m, a, a(Z, m)) and the observed frequency 
/ = (1 + 


The sources are stored into frequency bins of length 
A/ = 10 Hz, with central frequencies fi = lOi -I- 5 Hz 
for i=1...500. The energy density at the frequency fi is 
given by the discrete sum of the individual contribution 
of the Ni sources in the bin: 




SttG N f \ / rz\^^^BH{Z){l + z) 


(24) 

The spectrum for a = 0.2 and e = 10 ® is shown in Fig¬ 
ure [5] and compared to Models 1 and 2. The fluctuation 
at high frequencies is due to the small number sources 
in the simulation at such high frequencies. The above 
procedure is repeated for a grid of points in the e — a pa¬ 
rameter space, producing Hgw(/) for each point in this 
grid. The gravitational-wave spectra are then compared 
to the projected sensitivities of the second and third gen¬ 
eration detectors, similarly to Models 1 and 2, using the 
likelihood formalism (c.f. Equation [T5]) . 

Finally, to assess the importance of possible black hole 
progenitors of masses beyond our assumed upper limit of 
150Mq, we repeat the above calculations extending the 
IMF up to 300Mq and IOOOMq under the assumption 
that all stars above 300Mq produce IOOMq black holes. 

The 95% sensitivities of future detectors are shown in 
Figure [9l in the e — a plane. Note that the sensitivity 
contours are about 2x lower in e relative to Model 2, 
which is the result of the more careful treatment of the 
dependence of a on the metallicity, including the PISNe 
mass region. Changing the value of maximum progeni¬ 
tor mass again has a significant impact on detectability 
of this model, as already observed in Model 1. The two 
IMF models considered in this calculation agree to within 
a factor of about 2. We have verified that the two dif¬ 
ferent SFR models considered in this paper yield nearly 
identical results, as already observed in Models 1 and 2. 

Finally, as noted in the case of Model 1, the realistic 
values of e ^ 10“^ —10“® may be reached by Einstein 

Telescope, especially in the case if massive (> 300Mq) 
black hole progenitors exist, as suggested by recent ob¬ 
servations. The second generation detectors may pro be 
the more exotic scenarios associated with e ^ 0.01 1981. 


CONCLUSIONS 

In this paper we have studied the stochastic gravita¬ 
tional wave background generated by stellar core col¬ 
lapse events occurring throughout the universe. While 
the core collapse process is likely associated with several 
mechanisms for GW production, we focused on the most 
tractable one, namely the ringdown of the newly formed 
black hole, following the collapse itself. While this mech¬ 
anism has been studied in the past, our purpose here is 
to revisit it in more detail, and study the detectability 
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FIG. 9: 95% confidence sensitivity contours for the Advanced LIGO (black) and Einstein Telescope (gray) and for three values 
of mmax = I 5 OM 0 (solid), 300Mq (dashed), and IOOOM 0 (dot-dashed) are shown in the e — a plane for the Model 3, assuming 
the Salpeter IMF (left plot) and the Koupra IMF (right plot). 


of the corresponding background by the upcoming sec¬ 
ond and third generation gravitational-wave detectors. 
We considered three variations of the model. Our Model 
1 assumes that the dominant gravitational-wave mode 
during the ringdown is / = 2, and it integrates this sig¬ 
nal across all core collapse events that yield black holes. 
For this purpose, we used the most recent models of star 
formation rate that are consistent with the metallicity 
observations as well as with the CMB-based constraints 
on the reionization redshift and the optical depth param¬ 
eters. We scan the parameter space of the model and 
determine the part of the parameter space that is acces¬ 
sible to future detectors. In our Model 2, we attempt 
to include the effect of metallicity (which evolves with 
redshift) on the mass of the newly produced black hole. 
In this model we treat this effect in an average sense, 
allowing us to compute the gravitational-wave spectrum 
analytically. In Model 3, we go a step further and instead 
of an analytic calculation, we perform Monte Carlo sim¬ 
ulations, drawing the stars from the appropriate redshift 
and mass distributions determined by the Star Track sim¬ 
ulations. 

Remarkably, the three models agree with each other to 
within a factor of 2, indicating their robustness. How¬ 
ever, the required efficiency of gravitational-wave pro¬ 
duction by the ringdown is relatively high. In particu¬ 
lar, it is unlikely that the second generation detectors 
will observe this background since the required efficiency 
of gravitational-wave production is e ^ 10“"^ or higher, 
which is not supported by simulations. The third gen¬ 
eration detectors, however, would require e ^ 10~® or 
higher, which is more realistic. We emphasize, however, 
that we have considered only one of the several mecha¬ 
nisms for gravitational-wave production in core collapse 
processes. Future studies should attempt to repeat sim¬ 
ilar studies for other mechanisms—while other mecha¬ 
nisms are more difficult to model (due to the complexity 
of the relevant processes), they are also likely to produce 


a stronger gravitational-wave background. Furthermore, 
since most massive stars reside in binaries, we expect that 
inclusion of binary systems in this calculation may have a 
significant impact on the gravitational-wave background 
estimate. We plan to include this step in a follow-up 
study. 


Acknowledgments 

The work K.A.O. was supported in part by DOE 
grant de-sc0011842 at the University of Minnesota. The 
work of K.C. and V.M. was supported in part by 
NSF grant PHY-1204944 at the University of Min¬ 
nesota. KB and WG evolutionary modeling was par¬ 
tially sponsored by the NGN grant Sonata Bis 2 (DEC- 
2012/07/E/ST9/01360) and Polish Science Foundation 
(FNP) ’’Master 2013” Subsidy. The work of EV has been 
carried out at the ILP LABEX (under reference ANR- 
lO-LABX-63) supported by French state funds man¬ 
aged by the ANR within the Investissements dAvenir 
programme under reference ANR-ll-IDEX-0004-02, and 
was also sponsored by the French Agence Nationale pour 
la Recherche (ANR) via the grant VACOUL (ANR-2010- 
Blan-0510-01). Iigo-pl500077, 


[1] L. P. Grishchuk, Sov. Phys. JETP 40, 409 (1975). 

[2] A. A. Starobinskii, JETP Lett. 30, 682 (1979). 

[3] R. Eastlier and E. A. Lim, JGAP 0604, 010 (2006). 

[4] N. Barnaby, E. Pajer, and M. Peloso, Phys. Rev. D 85, 
023525 (2012). 

[5] R. R. Galdwell and B. Allen, Phys. Rev. D 45, 3447 
(1992). 

[6] T. Damour and A. Vilenkin, Phys. Rev. Lett. 85, 3761 

( 2000 ). 

[7] T. Damour and A. Vilenkin, Phys. Rev. D 71, 063510 
(2005). 















10 


[8] X. Siemens, V. Mandic, and J. Creighton, Phys. Rev. 
Lett. 98, 111101 (2007). 

[9] S. Olmez, V. Mandic, and X. Siemens, Phys. Rev. D 81, 
104028 (2010). 

[10] S. Olmez, V. Mandic, and X. Siemens, J. Cosm. Astrop. 
Phys. 07 , 009 (2012). 

[11] M. Gasperini and G. Veneziano, Astropart. Phys. 1 , 317 
(1993). 

[12] V. Mandic and A. Buonanno, Phys. Rev. D 73 , 063008 
(2006). 

[13] T. Regimbau, Res. Astr. Astrop. 11 , 369 (2011). 

[14] E. S. Phinney, Astrophys. J. Lett. 380, L17 (1991). 

[15] D. 1. Kosenko and K. A. Postnov, Astron. & Astrop. 
336 , 786 (1998). 

[16] T. Regimbau and J. A. de Freitas Pacheco, Astrophys. 
J. 642, 455 (2006), gr-qc/0512008. 

[17] X.-J. Zhu, E. Howell, T. Regimbau, D. Blair, and Z.-H. 
Zhu, Astrophys. J. 739 , 86 (2011). 

[18] P. A. Rosado, Phys. Rev. D 84, 084004 (2011), 
1106.5795. 

[19] S. Marassi, R. Schneider, G. Corvino, V. Ferrari, and 
S. P. Zwart, Phys. Rev. D 84, 124037 (2011), 1111.6125. 

[20] C. Wu, V. Mandic, and T. Regimbau, Phys. Rev. D 85, 
104024 (2012). 

[21] X.-J. Zhu, E. J. Howell, D. G. Blair, and Z.-H. Zhu, 
Monthly Notices of the Royal astronomical Society 431 , 
882 (2013), 1209.0595. 

[22] E. F. D. Evangelista and J. G. N. de Araujo, Modern 
Physics Letters A 28, 1350174 (2013), 1504.04300. 

[23] 1. Kowalska-Leszczynska, T. Regimbau, T. Bulik, 
M. Dominik, and K. Belczynski, A&A 574, A58 (2015), 
1205.4621. 

[24] E. F. D. Evangelista and J. C. N. de Araujo, Mon. Not. 
Roy. Astron. Soc. 449 , 2700 (2015), 1504.02700. 

[25] T. Regimbau and J. A. de Freitas Pacheco, Astron. and 
Astrophys. 376 , 381 (2001). 

[26] B. J. Owen, L. Lindblom, C. Gutler, B. F. Schutz, 
A. Vecchio, and N. Andersson, Phys. Rev. D 58 , 084020 
(1998). 

[27] V. Ferrari, S. Matarrese, and R. Schneider, Monthly No¬ 
tices of the Royal astronomical Society 303 , 258 (1999), 
astro-ph/9806357. 

[28] S. Chandrasekhar, ’’Ellipsoidal Figures of Equilibrium” 
(New Haven, Yale Univ. Press, 1969). 

[29] J. L. Houser, J. M. Centrella, and S. C. Smith, Phys. 
Rev. Lett. 72 , 1314 (1994). 

[30] D. Lai and S. L. Shapiro, Astrophys. J. 442 , 259 (1995). 

[31] E. Howell, D. Coward, R. Burman, D. Blair, and 
J. Gilmore, Monthly Notices of the Royal astronomical 
Society 351 , 1237 (2004). 

[32] X.-J. Zhu, X.-L. Fan, and Z.-H. Zhu, Astrophys. J. 729, 
59 (2011), 1102.2786. 

[33] P. A. Rosado, Phys. Rev. D 86, 104007 (2012), 
1206.1330. 

[34] P. D. Lasky, M. F. Bennett, and A. Melatos, Phys. Rev. 
D 87 , 063004 (2013), 1302.6033. 

[35] G. Gutler, Phys. Rev. D 66, 084025 (2002). 

[36] T. Regimbau and J. A. de Freitas Pacheco, Astronomy 
and Astrophysics 447 , 1 (2006), astro-ph/0509880. 

[37] T. Regimbau and V. Mandic, Glass. Quant. Grav. 25 , 
184018 (2008). 

[38] S. Marassi, R. Giolfi, R. Schneider, L. Stella, and V. Fer¬ 
rari, Monthly Notices of the Royal astronomical Society 
411 , 2549 (2011), 1009.1240. 


[39] E. Howell, T. Regimbau, A. Gorsi, D. Goward, and 
R. Burman, Monthly Notices of the Royal astronomi¬ 
cal Society 410 , 2123 (2011), 1008.3941. 

[40] C.-J. Wu, V. Mandic, and T. Regimbau, Phys. Rev. D 
87, 042002 (2013). 

[41] P. Sandick, K. A. Olive, F. Daigne, and E. Vangioni, 
Phys. Rev. D 73, 104024 (2006). 

[42] A. J. Farmer and E. Phinney, Mon. Not. Roy. Astron. 
Soc. 346, 1197 (2003). 

[43] B. Abbott et ah, Phys. Rev. Lett. 95, 221101 (2005). 

[44] B. Abbott et ah, Astrophys. J. 659 , 918 (2007). 

[45] B. Abbott et ah. Nature 460, 990 (2009). 

[46] J. Aasi et ah, Phys. Rev. Lett. 113, 231101 (2014). 

[47] J. Aasi et ah, Phys. Rev. D 91, 022003 (2015). 

[48] B. Abbott et ah, Phys. Rev. D 76, 082003 (2007). 

[49] B. Abbott et ah, Phys. Rev. Lett. 107, 271102 (2011). 

[50] B. Abbott et ah, Nucl. Instr. Meth. A 517, 154 (2004). 

[51] B. Abbott et ah. Rep. Prog. Phys. 72, 076901 (2009). 

[52] F. Acernese et ah. Class. Quant. Grav. 23, S63 (2006). 

[53] V. Mandic, E. Thrane, S. Giampanis, and T. Regimbau, 
Phys. Rev. Lett. 109 , 171102 (2012). 

[54] S. G. Crowder, R. Namba, V. Mandic, S. Mukhoyama, 
and M. Peloso, Phys. Lett. B 726 , 66 (2013). 

[55] G. M. Harry (for the LIGO Scientific Gollaboration), 
Class. Quant. Grav. 27, 084006 (2010). 

[56] https : //wwwcascina. virgo . inf n. it/advirgo/docs . html 

[57] B. Willke et ah. Glass. Quant. Grav. 23, S207 (2006). 

[58] G. Affeldt et ah. Glass. Quant. Grav. 31, 224002 (2014). 

[59] Y. Aso and Y. Michimura and K. Somiya and M. Ando 
and O. Miyakawa and T. Sekiguchi and D. Tatsumi and 
and H. Yamamoto, Phys. Rev. D 88, 043007 (2013). 

[60] K. S. (for the KAGRA collaboration). Class. Quant. 
Grav. 29 , 124007 (2012). 

[61] T. E. science team, https://tds.ego-gw.it/ql/?c=7954 

( 2011 ). 

[62] C. Ott, G. Reisswig, E. Schnetter, E. O’Gonnor, 

U. Sperhake, F. Loffler, P. Diener, E. Abdikamalov, 

1. Hawke, and A. Burrows, Phys. Rev. Lett. 106, 161103 

( 2011 ). 

[63] C. Ott et ah, Astrophys. J. 768, 115 (2013). 

[64] B. Muller et ah, Astrophys. J. 768, 115 (2013). 

[65] V. Ferrari, S. Matarrese, and R. Schneider, Mon. Not. 
Roy. Astron. Soc. 303 , 247 (1999), astro-ph/9804259. 

[66] J. C. N. de Araujo, O. D. Miranda, and O. D. Aguiar, 
Phys. Rev. D 61 , 124015 (2000), astro-ph/0004395. 

[67] J. G. N. de Araujo, O. D. Miranda, and O. D. Aguiar, 
Mon. Not. Roy. Astron. Soc. 330 , 651 (2002), astro- 
ph/0202037. 

[68] J. de Araujo, O. Miranda, and O. Aguiar, Glass. Quant. 
Grav. 19 , 1335 (2002). 

[69] J. G. N. de Araujo, O. D. Miranda, and O. D. Aguiar, 
mnras 348 , 1373 (2004). 

[70] J. C. N. de Araujo, O. D. Miranda, and O. D. Aguiar, 
Classical and Quantum Gravity 21, 545 (2004). 

[71] A. Buonanno, G. Sigl, G. Raffelt, H. Janka, and 
E. Muller, Phys. Rev. D 72, 084001 (2005). 

[72] D. Coward, R. Burman, and D. Blair, Mon. Not. Roy. 
Astron. Soc. 329 , 411 (2002). 

[73] X.-J. Zhu, E. Howell, and D. Blair, Monthly Notices 
of the Royal Astronomical Society 409 , L132 (2010), 
1008.0472. 

[74] S. Marassi, R. Schneider, and V. Ferrari, Mon. Not. Roy. 
Astron. Soc. 398, 293 (2009). 

[75] F. Pacucci, A. Ferrara, and S. Marassi, 


11 


arXiv: 1502.04125 (2015). 

[76] F. Echeverria, Phys. Rev. D 40, 3194 (1989). 

[77] R. Stark and T. Piran, Phys. Rev. Lett. 55, 891 (1985). 

[78] K. Belczynski et al., http: //www. syntheticuniverse . org 

[79] S. Lilly, O. L. Fevre, F. Hammer, and D. Crampton, 
Astrop. J. 460, LI (1996). 

[80] A. M. Hopkins and J. Beacom, Astrop. J. 651, 142 
(2006). 

[81] M. A. Fardal, N. Katz, D. H. Weinberg, and R. Dav, 
Mon. Not. Roy. Astron. Soc. 379 , 985 (2007). 

[82] S. M. Wilkins, N. Trentham, and A. Hopkins, 
arXiv:0803.4024 (2008). 

[83] K. Nagamine, J. P. ©striker, M. Fukugita, and R. Cen, 
Astrop. J. 653 , 881 (2006). 

[84] L. Hernquist and V. Springel, Mon. Not. Roy. Astron. 
Soc. 341 , 1253 (2003). 

[85] P. Behroozi, R. Wechsler, and C. Conroy, Astrophys. J. 
770 , 57 (2013). 

[86] P. Oesch et al., Astrophys. J. 786, 108 (2014). 

[87] P. Oesch et ah, arXiv: 1409.1228 (2014). 

[88] B. Robertson and R. Ellis, Astrophys. J. 744, 95 (2012). 

[89] F. Y. Wang, A&A 556, A90 (2013). 

[90] M. Kistler, H. Yuksel, and A. Hopkins, arXiv:1305.1630 
(2013). 

[91] M. Trenti, R. Perna, and S. Tacchella, Astrophys. J. 
Lett. 773 , 22 (2013). 

[92] P. Behroozi and J. Silk, Astrophys. J. 799, 32 (2015). 

[93] E. Vangioni, K. Olive, T. Prestegard, J. Silk, P. Petit- 
jean, and V. Mandic, arXiv:1409.2462 (2014). 

[94] G. Hinshaw et ah, Astrophys. J. Supl 208, 19 (2013). 

[95] P. A. Crowther, O. Schnurr, R. Hirschi, N. Yusof, R. J. 
Parker, S. P. Goodwin, and H. A. Kassim, Mon. Not. 
Roy. Astron. Soc. 408 , 731 (2010), 1007.3284. 

[96] https : //dec. ligo . org/LIG0-T0900288/public 

[97] K. Belczynski, T. Bulik, C. Fryer, A. Ruiter, F. Valsec- 
chi, J. Vink, and J. Hurley, Astrophys. J. 714, 1217 
( 2010 ). 

[98] S. Kobayashi and P. Meszaros, Astrophys. J. 589 , 861 
(2003). 

[99] L. Baiotti and L. Rezzolla, Phys. Rev. Lett. 97, 141101 
(2006). 


[100] S. Woosley and T. Weaver, Astrophys. J. Supl 101, 181 
(1995). 

[101] K. Belczynski, V. Kalogera, and T. Bulik, Astrophys. J. 
572, 407 (2002), arXiv:astro-ph/0111452. 

[102] K. Belczynski, V. Kalogera, F. A. Rasio, R. E. Taam, 
A. Zezas, T. Bulik, T. J. Maccarone, and N. Ivanova, 
Astrophys. J. Supl 174, 223 (2008), arXiv:astro- 
ph/0511811. 

[103] J. R. Hurley, O. R. Pols, and G. A. Tout, Mon. Not. Roy. 
Astron. Soc. 315, 543 (2000), arXiv:astro-ph/0001295. 

[104] M. Dominik, K. Belczynski, G. Fryer, D. E. Holz, 
E. Berti, T. Bulik, 1. Mandel, and R. O’Shaughnessy, 
Astrophys. J. 759, 52 (2012), 1202.4901. 

[105] H. Sana, S. E. de Mink, A. de Koter, N. Langer, C. J. 
Evans, M. Gieles, E. Gosset, R. G. Izzard, J.-B. Le 
Bouquin, and E. R. N. Schneider, Science 337, 444 
(2012), 1207.6397. 

[106] J. S. Vink, A. de Koter, and H. J. G. L. M. Earners, 
A&A 369, 574 (2001), astro-ph/0101509. 

[107] G. L. Fryer, K. Belczynski, G. Wiktorowicz, M. Do¬ 
minik, V. Kalogera, and D. E. Holz, Astrophys. J. 749, 
91 (2012), 1110.1726. 

[108] K. Belczynski, G. Wiktorowicz, C. L. Fryer, D. E. 
Holz, and V. Kalogera, Astrophys. J. 757, 91 (2012), 
1110.1635. 

[109] F. Ozel, D. Psaltis, R. Narayan, and J. E. McGlintock, 
Astrophys. J. 725, 1918 (2010), 1006.2834. 

[110] G. D. Bailyn, R. K. Jain, P. Coppi, and J. A. Orosz, 
Astrophys. J. 499, 367 (1998), astro-ph/9708032. 

[111] G. L. Fryer, S. E. Woosley, and A. Heger, api 550, 372 
(2001), astro-ph/0007176. 

[112] N. Yusof, R. Hirschi, G. Meynet, P. A. Crowther, S. Ek- 
strom, U. Frischknecht, C. Georgy, H. Abu Kassim, and 
O. Schnurr, Mon. Not. Roy. Astron. Soc. 433, 1114 
(2013), 1305.2099. 

[113] K. Belczynski, A. Buonanno, M. Cantiello, G. L. Fryer, 
D. E. Holz, 1. Mandel, M. G. Miller, and M. Walczak, 
Astrophys. J. 789, 120 (2014), 1403.0677. 

[114] P. Kroupa and G. Weidner, Astrophys. J. 598, 1076 
(2003), astro-ph/0308356. 


